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Abstract 

Recently, out-of-equiUbrium field theory has been studied using approxima- 
tions based on truncations of the 2PI effective action. Although results are 
promising, the convergence of subsequent orders of the approximation is diffi- 
cult to get a handle on, mainly because, generically, only the lowest non-trivial 
order is currently numerically tractable. We study one specific case, the broken 
phase of the (j)'^ model, where the existence of an effective three point vertex 
makes it possible to compare a lowest and next-to-lowest non-trivial order. 

1 Introduction 

Out-of-equilibrium quantum fields can be treated using real-time equations of motion based on trun- 
cations of the 2PI effective action fT "21 (based on the formalism of jSj |4l). The equations are writ- 
ten in terms of the (connected) propagator (0(x, t)(?!)(0, t))c = G{x,t,t') = F(x, t, t') — sign(t — 
t')ip{yL,t,t')/2, the mean field {4>iO,t)) = (p{t) and the self-energies J:{x,t,t') = Sioc(0, t)(5(t - 
t')5{x) + S^(x, t, t') - sign(t - t')iT.plx, t, t'), and read 



dt" / (i^zS''(z,t,f")F(x-z,t",t' 



(1) 



^* dt" J d^z S^(z, t, t")p{x - z, t", t'), 



dt" / d^zJ:P{z,t,t")p{x-z,t",t'), 



(5^[(/>, G] 



(2) 



with M|?jj(t) = m? + Sioc(0, t). A ^-derivable approximation results from specifying a truncation of 
the infinite set of 2PI diagrams $ contributing to the self energy through. 



S(x,t,0 



5G{x,t',t) 

Initial conditions are set by specifying an initial gaussian density matrix through 



(3) 



({0o(k),</)o(-k)}) = (nfe + l/2)M, ({7ro(k),^o(-k)}) = K + 1/2K, ({7ro(k), 0o(-k})) = 0, 

(4) 

for some choice of distribution function and dispersion relation cok, and implementing the relations 
([0o(k),<^o(-k)]) = 0, ([^o(k),0o(-k)]) = 1. 

Two expansions are on the market; one in the 0{N) model in powers of f5"2"6"7l which has 
now become commonly used, and one in the coupling LL<8.^> which we will use here. Although it may 
potentially be less broadly applicable, it is straightforward to use and easy to understand. It is important 
to stress that the equations are self-consistent in a similar way to a gap equations, and that they resum 
infinite sets of diagrams. 
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Fig. 2: The included diagrams to 0(A ) in the effective action (first line). Lines are full propagators G, small circles the mean 
field (j). In the symmetric phase (0 — 0) the first, second and fourth diagram vanish. In the broken phase </> oc A~^/^, they are 
all there, as well as two more at 0{\^) (in parenthesis), which we do not include . Below, the derived self-energy contributions 
for the mean field (second line) and propagator (third line). 



2 This work 

We study truncations of the coupling expansion of (f>'^ theory, with the action 

The model has two phases: the "symmetric" {m? > 0), where ((/>) = 0, and the "broken" (m^ < 0), 
with {(/))= V = a/6 m'-^/X. Figure 1 shows the diagrams at the level of the action and the resulting 
contributions to the self-energies in orders of A in either phase. The expansion is truncated at the order 
which is currently numerically tractable; including further vertices leads to additional nested space-time 
"memory integrals" on the right hand side of the equations of motion At O(A^), two additional 
diagrams exist for non-zero mean field. They are also not numerically tractable, and we consider the one 
(Basketball) diagram included here to represent the effect of the full O(A^) contribution. The aim is to 
determine whether one order dominates the next. 

In the symmetric phase, only the "leaf" diagram contributes at 0(A), leading to the Hartree ap- 
proximation. We will consider this a trivial approximation, since it is simply a (time-dependent) con- 
tribution to the effective mass (ReS = Sioc(0,t)), and it is included in the quantity Af^jj(t). Mode 
damping, equilibration and thermalization, processes of interest in out-of-equilibrium field theory, are 
the result of a non-zero imaginary part of the self-energy and in the symmetric phase, this appears at 
O(A^), in the last (Basketball) diagram. Detailed study of thermalization and equilibration in this trun- 
cation in 1, 2 and 3-1-1 dimensions can be found in JQJ. In particular, the damping rate can be 
calculated in perturbation theory, the result being qualitatively compatible with the 2PI resummed result 
lITOl . Apparently, the resummation has only a small, quantitative effect in this case. 

In the broken phase, the mean field leads to an effective three point vertex, and because {(/)) is 
formally 0(A^^/^), an additional 2-loop diagram appears at order 0(A) (Two-loop) in the effective 
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Fig. 3: Distribution function 71^ + 1/2 at various times. Around time mt — 500, tlie curve is smootli and approaciiing a 
Bose-Einstein distribution. 



action. Although in perturbation theory it does not lead to on-shell damping fTO\, resummation from the 
self-consistent evolution equations leads to a non-zero ImS. In this case, we are able to compare two 
successive truncations of ImS and estimate the quality of the expansion. 

3 Results 

Figure 2 shows the evolution of the effective particle number + 1/2 = y^F(k^T7F)^^rF(k7M0t=t' 
in time. We start with a "top-hat" (Tl) initial condition which is very out-of-equilibrium, and we can 
see a progression to a smooth Bose-Einstein-like distribution at time mt ~ 500. The dark (black) dots 
are O(A^), the light (green) 0(A). It is clear, that not only does the resummation lead to a quantitative 
effect (the non-zero ImS at 0(A)), but the time scales of equilibration are also very similar. Here, A = 1. 
This constitutes kinetic equilibration, and will (at small coupling) lead to a Bose-Einstein distribution 
with an effective temperature and chemical potential. The (p^ model has no conserved charges, and the 
asymptotic equilibrium state has zero chemical potential. To reach it requires chemical equilibration to 
take place as well. 

Figure 3 shows ln(l + l/n^) vs. Uk = y/dtdt'F{k, t, t') /i='(k, t, t')t=t/ at mt = 1000. A Bose- 
Einstein distribution is then a straight line with slope T^^ and intercept fich/T. Overlaid are the results 
for Bose-Einstein distributions with zero chemical potential. Clearly, the low-momentum range has 
equilibrated kinetically, but not (yet) chemically. It is well known that at O(A^), chemical equilibration 
takes place (see for instance im i^ lTni ). In figure 4 we plot the total particle number A'tot = J d^^n-k 
normalized to the initial value. In both approximations, there is a change of particle number. Although 
the behaviour is not the same, presumably they will both lead to chemical equiUbration. 
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Fig. 4: The distribution function ln(l + l/rife) at time mt — 1000. The low momentum range has kinetically equilibrated. The 
straight line are Bose-Einstein distributions at various temperatures and zero chemical potential. 

4 Conclusions 

Successive, non-trivial truncations of the coupling expansion of the 2PI/$-derivable approximation have 
been compared. In the broken phase of (j)'^ theory, kinetic and chemical equilibration is included at 0(A), 
at a rate very close to what we find at 0(A) + O(A^). At this moderate coupling, the expansion seems 
to be well behaved. In addition, we found that the relevant diagram at 0(A), which perturbatively does 
not lead to on-shell damping, becomes the dominant contribution when resummed by the self-consistent 
evolution equations. 

Further investigation of the convergence of expansions of the 2PI effective action is needed. Be- 
cause of the numerical effort required, it may be necessary to consider simpler systems, such as quantum 
mechanics, as a testing ground. For further details and a study of equilibration in the symmetric phase, 
see flOl . 
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Fig. 5: The total particle number in time. 
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